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We investigate the dynamics of a fluid of dipolar hard 
spheres in its liquid and glassy phase, with emphasis on the 
microscopic time or frequency regime. This system shows 
rather different glass transition scenarios related to its rich 
equilibrium behavior which ranges from a simple hard sphere 
fluid to a long range ferroelectric orientational order. In the 
liquid phase close to the ideal glass transition line and in 
the glassy regime a medium range orientational order occurs 
leading to a softening of an orientational mode. To investigate 
the role of this mode we use the molecular mode-coupling 
equations to calculate the spectra (ji'i'milt ^) ^i^d Xi'milj In 
the center of mass spectra 000 (Qi ^^'^ Xooil,^) we found 
besides a high frequency peak at uj^f a peak at uiop, about 
one decade below Whf- oJop has almost no q-dependence and 
exhibits an "isotope" effect u)op oc , with 1 the moment 

of inertia. We give evidence that the existence of this peak is 
related to the occurrence of the medium ranged orientational 
order. It is shown that some of these feature also exist for 
schematic mode coupling models. 

61.25.Em, 64.70.Pf, 61.20.Lc 



I. INTRODUCTION 

The dynamical properties of liquids and glasses still 
are a challenging problem. In the vicinity of the glass 
transition the frequency (or time) range can be decom- 
posed into at least three different domains, the a, /3 and 
microscopic regime. The first one describes the struc- 
tural relaxation which dramatically slows down when the 
glass transition is approached from above. In the ideal- 
ized mode-coupling theory (MCT) for simple liquids 0-^] 
and for molecular systems 1^-0] it even stops at a critical 
temperature T^- Probably the most interesting result of 
MCT is the existence of the so-called /3-relaxation which 
describes the dynamics within a cage of particles above 
and below Tc. The corresponding /?- frequency scale is 
much larger than that for the a-relaxation. At still higher 
frequencies there are vibrational and librational motion 
which constitute the microscopic regime. 

One may say that most of the attention in the field of 
glassy dynamics during the last fifteen years has been de- 
voted to the a- and /3 -relaxation. This activity has been 



mainly stimulated by MCT which has predicted in these 
two regimes scaling laws with a diverging a- and /9-time 
scale. These predictions have been tested intensively by 
experiments and numerical simulations. A satisfactory 
agreement has been found for many glass forming systems 
P-jrTl,P|jl2|Jl^j7|,p^ . Experimental and simulational re- 
sults do not exhibit any singular or crossover behavior 
for microscopic frequencies lo, i.e. for lo>\ THZ. Never- 
theless in that regime an interesting phenomenon occurs 
in most glasses but not in crystals and colloidal glasses, 
which is the so called boson-peak. Indications for this 
peak came from two different sides. First, the low tem- 
perature specific heat c{T)/cd{T) scaled by the phonon 
contribution cd{T) oc shows apeak at about 10 K 
for several glass formers (see e.g. ||l^). This excess with 
respect to cd{T) around 10 K implies the existence of 
additional excitations besides the long wavelength acous- 
tic phonons. Second, Raman spectra I{uj) compared to 
the phonon contribution Id{<^) oc exhibit at about 1 
THZ an excess, as well (see e.g. |l^). Since the tem- 
perature dependence of the excess intensity scales with 
the Bose-distribution fimction nsiT), the peak is called 
boson-peak. 

That these two observations might have a common ori- 
gin was first shown by Buchenau et al. |^,|9|. For vit- 
reous silica these authors also found an excess with re- 
spect to goi^) oc for the vibrational density of states 
g{Lo) determined from inelastic neutron scattering (INS) 
data. Using g{uj) to calculate c{T) / cd{T) led to a good 
agreement with the result from heat capacity measure- 
ments. Particularly, the peak positions for both results 
coincided. Although the boson-peak does not seem to 
posses any singular ld- or T-dependence it is a universal 
phenomenon for all systems with exclusion of glass form- 
ing colloids, in the sense that it appears more or less for 
almost all glass formers. For, e.g. LiCl solutions jl|] and 
Orthoterphenyl , it has been stressed that the boson- 
peak, together with the narrowing of a central peak at 
w = 0, develops continously from the liquid to the glassy 
phase. 

Despite considerable experimental and numerical ef- 
forts its microscopic origin is still not satisfactorily under- 
stood. Two reasons might be responsible for that. First, 
e.g. in case of light scattering the precise connection 
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between the measured quantity and the basic theoreti- 
cal objects, the time- or frequency-dependent site-site or 
molecular correlation functions for molecular liquids, is 
not known. For instance, it has been shown that several 
coupling mechanism between light and distinct modes of 
liquid ZnCl2 exist which have different w-dependent cou- 
pling constants | |2l| ] . This may complicate the determina- 
tion of g{w) from Raman spectra. Second, the relation- 
ship between the various t- or u- dependent correlators 
and microscopic modes obtained from a diagonalization 
of the dynamical matrix is not obvious. In addition, these 
correlators can only be calculated analytically under se- 
rious approximations. Apparently, numerical investiga- 
tions represent a powerful tool since their microscopic 
nature allows to calculate both, the correlators or spec- 
tra and under certain conditions the microscopic modes. 
The possible character of these modes range between: 



vibrational 
extended 
propagating 
harmonic 
longitudinal 
coherent 
translational 
acoustic 



relaxational 

localized 

non-propagating 

anharmonic 

transversal 

random 

rotational 

optic 



Since these feature can occur in combinations, the com- 
plexity of the problem becomes obvious. In addition, 
different experimental approaches to a specific material 
or a specific measurement of different type of glass for- 
mers, e.g. strong or fragile ones, may exhibit different 
features of the same phenomenon. 

Without demanding completeness let us shortly review 
the present status. Inelastic X-ray scattering (IXS) on 
v-Si02 gave evidence that propagating acoustic sound 
waves exist even above to bp-, the position of the boson- 
peak . This has let these authors to conclude that the 
propagating modes are also involved in the boson-peak 
itself. On the other side a crossover at lobp from prop- 
agating to localized (strongly scattered) acoustic modes 
was deduced for v-Si02 from INS and IX S P3[. Com- 
ments on this controversy are given in Ref . ||24|-p6[ . How- 
ever, interpretation of IXS-experiments which are free of 
any model have recently strengthened at least the fact 
that there are propagating modes above lubp p7| , p8t . 
MD-simulations for v-Si02 have given evidence that the 
boson-peak modes can not strongly be localized and that 
there is a contribution of transverse propagating modes 
in the boson-peak regime j2^,^^. IXS for LiCl-solutions 
and glycerol seem to detect propagating, longitudinal 
modes at ljbp ||3^, whereas a MD-simulation for II2O 
shows a mixing of propagating longitudinal and transver- 
sal modes, at least for q large enough {q > Anm^^) |32|. 

Further important information about the nature of the 
boson-peak modes comes from a normal mode analysis 



for a system of soft spheres [Q and for Si02 deep in its 
glass phase |3^-|37|. There it has been found that these 
modes are harmonic and quasi localized p3| , ^ . They 
occur due to hybridization of localized low frequency op- 
tic modes with propagating acoustic states [3^ ] and are 
extended but non-propagating. In addition each normal 
mode has a coherent and a random component. The lat- 
ter finding is consistent with other numerical results for 
a Lennard- Jones liquid |§, liquid ZnCl2 and 
conclusions drawn from experimental data for various 
glass formers Since ZnCl2 has been considered 

in its liquid phase an instantaneous normal mode analy- 
sis is feasible, only. Restriction to the stable modes accu- 
rately reproduces in the microscopic frequency regime the 
dynamical structure factor S{q, uj) obtained from a MD- 
simulation pO| . This implies a harmonic mode charac- 
ter at temperatures even above the MCT glass-transition 
temperature Tc- 

There is not much analytical work on the boson peak. 
For a harmonic crystal with random spring constants it 
has been shown that an excess density of states 

follows. It has recently been stressed that the ap- 
proach in has some shortcomings which may be re- 
moved by use of off-lattice models Q. Another theoret- 
ical framework is the soft potential model p7|-p0t which 
allows to describe the low temperature anomalies below 
IK and those at 10-20 K, which includes the excess den- 
sity of states. 

As already mentioned above, MCT has been rather 
successful in describing the a- and /3-dynamics. 
Schematic MCT-models [0 where the wave number de- 
pendence has been neglected were also used to describe 
experimental spectra in the microscopic regime, includ- 
ing the boson peak | |Tl| , ^l| -^6t . Whereas the results in 
Ref. [|n| (cf. Fig. 6) and Ref. || (cf. Fig. 4) yield 
spectra which are reminiscent to boson peak spectra the 
line shape does not come out, satisfactorily. However, 
recently a detailed MCT-investigation for a glass of hard 
spheres was performed, including the q- dependence \ ^t\ . 
Besides a high frequency peak at lohf an additional peak 
at ujAOP, about one decade below luhf, has been seen for 
a volume fraction ip = 0.6. This peak, which strongly re- 
sembles that in Refs. (Fig. 6) and [|| (Fig. @), origi- 
nates from the distribution of harmonic oscillators within 
the cages and has been called anomalous oscillation peak 



(AOP). As explicitely demonstrated |57| it shares many 
features with the boson peak. 

In contrast to a system of hard spheres or binary van 
der Waals liquids, molecular liquids also have orienta- 
tional degrees of freedom. The question which one may 
ask is: Does the boson peak also involve orientational 
motion? Indeed one of its first interpretations for v-Si02 
were coupled rotations of Si04-tetrahedra jl^. Such 
an interpretation has been supported by MD-simulations 
for v-Si02 p5| , |58[ | and ZnCl2 |Q and by neutron scat- 
tering experiments, proving non-sound-like contributions 
around ujbp ]59[ |. Dielectric loss measurements probe 
the orientational dynamics, only. Since these measure- 
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ments, e.g. for Glycerol and Propylene Carbonate, also 
exhibit a boson peak |6^j6l[] gives additional evidence 
that this peak may also be related to the orientational 
degrees of freedom. The role of orientational modes be- 
comes even more clear from the experiments on ethanol 
1^,^. Around 100 K, Ethanol can occur in several 
phases: glass phase, orientational glass phase, crystalline 
phase and a rotator phase. The center of mass posi- 
tions of the molecules are frozen in an amorphous struc- 
ture for the glass phase, and in a crystalline structure 
for the three other phases. The orientational dynamic is 
non-ergodic for the glass and the orientational glass and 
ergodic for the crystal and the rotator phase. INS has 
shown the existence of a boson peak for the orientational 
glass which does not differ much from that in the glass. 
Even in the rotator phase there is a boson peak, how- 
ever, shifted to lower frequencies. Similarly to structural 
glasses the orientational glass phase also exhibits an ex- 
cess in the specific heat c{T)/co{T) |62|. These findings 
suggest that these excess modes are primarily related to 
the orientational degrees of freedom. 

In the present paper we will not calculate the density 
of states, but the susceptibility and correlation spectra 
of the collective dynamics. These quantities show the 
appearance of an extra peak about a decade below the 
high frequency peak, where the former also originates 
from the orientational degrees of motion. 

The outline of the paper is as follows. In Sec. ||, we 
will shortly review the derivation of the equation of mo- 
tion for the relevant correlators of rigid, linear molecules 
and the mode-coupling approximation. In addition, we 
will discuss the linearized equations. The molecular 
mode-coupling equations are solved numerically for a liq- 
uid of dipolar hard spheres in Sec. [IE. It is shown in 



Sec. Ill that the qualitative features of the susceptibility 



spectra obtained for dipolar hard spheres can be derived 
from a schematic model. Finally, Section ^ contains a 
discussion of the results and some conclusions. 



II. EQUATIONS OF MOTION 

In the first part of this section we will present the equa- 
tions of motion for the most relevant correlation functions 
of a molecular liquid. The second part contains a dis- 
cussion of the corresponding linearized equations, which 
yields information on the microscopic time scale and in 
the third part we will shortly review the mode-coupling 
approximation. 



A. Mori equation 

We restrict ourselves to a system of N rigid and lin- 
ear molecules with mass M and moment of inertia /. 
There are two possibilities to describe molecular liquids: 
a. site- site or a moZecMZar representation p4| . The latter. 



which will be chosen here, decomposes the 5N degrees 
of freedom into 3A^ translational and 2N orientational 
ones. The starting point is the microscopic tensorial one- 
particle density mode 



N 



Pli 



Mt) = ^e^?--^W Yiminnit)) (1) 



and the corresponding translational and rotational cur- 
rent density mode 



TV 



n=l 

for a = T and R, respectively. Xn{t) and f2„(t) = 
{9n{t), ipn{t)) are the center of mass position and the ori- 
entation (specified by the two angles On and ipn) of the 
n-th molecule at time t, respectively. Yim{0,ip) are the 
spherical harmonics with / = 0, 1, 2, —l<m<l and 
the velocities are defined by 



Xn{t) , a = T 
ujn{t) , a = i? 



(3) 



where is the corresponding angular velocity, pim 

and j;"^ are related to each other by the continuity equa- 
tion: 



pim {q, i) = i^ q? iq) jim {q, t) 



(4) 



a=T.,R 



with: 



jLM t) = ^q- jfnM t), q^\q] (5a) 
jL {q, t) = ^T^f^^ • ^^^r. {q. t) (5b) 

the "longitudinal" translational and rotational current 
density and 



qriq) 



q , a — T, \fl 



(6) 



L denotes the operator of angular momentum. 

The time- and {q; Im, Z'm')-dependent correlation func- 
tions: 

Slra, I'm' (?, = ^ Plm (?) t)5pVm' {q, 0)) (7) 

and 



J^nli'm' (q, t)^j^ {jZ iq, t) 3vm' {q, o)) (8) 

are of particular experimental and theoretical impor- 
tance, at least for I = I' < 2. Here, Spim{q,t) = 
Pimiq,t) — {pim{q,t)) is the corresponding fluctuation of 
the density. S{q,t) = {Sim,i'm'iq,t)) form a complete 
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set for any time dependent two-point density correla- 
tor. Particularly, the partial dynamical structure fac- 
tors in a site-site description are linear superpositions of 
the molecular correlators Sim,i'm'i<l, t), but not vice versa 
dssf . In addition, we introduce the density-current den- 
sity correlator: 



Jlr. 



X, y, z as slow variables. This has been done 



(9) 



Mori-Zwanzig formalism has been used to derive equa- 
tions of motion for S{q,t) for a single linear molecule in 
a liquid of isotropic particles |^ and for a molecular liq- 
uid of linear ^ and arbitrary molecules Hj. Similar work 
has been done in a site-site description pq] . A compari- 
son between the tensorial and the site-site mode coupling 
theory has recently been performed for a single dumbell 
in a liquid of hard spheres ||67|| . 

Choosing Spim and both scalar current densities 
Jim , Oi = T,R as slow variables, the Mori-Zwanzig for- 
malism leads to the following set of equations: 

S(g,t) -f i^q"S"((f,t) = (10a) 

a 

S"iq, t) + zq"J"(q)S-l (g, 0) S(q, t) + 
t 

J" {q) J dt' m""' {q, t-t') S"' (g, t')^Q (10b) 

a' 

with initial conditions: 



S(g,0) 
S(g,0) 



0= 



(9,0) 



S"(Q,0) = ~*J(q)"q" 



and 



(11) 

(12) 
(13) 

(14) 



(J" (9));™.'™' ^ JZj-m-{<^ = -j^5^^,6u'5ram' (15) 

Here, S{q) and J"" [q] denote, respectively, the static 
density and current density correlation matrix and 



M, a = T 
I, a = R 



(16) 



Note first that these equations are completely equiv- 
alent to the set of equations in [^. However, they are 
given in a different representation, which makes its nu- 
merical solution convenient. Second, this set of equations 
is still exact, but needs an expression for the memory ker- 
nels rnl^i,^, {q,t). This is where approximations come 
in. Their nature depends strongly on the physical sit- 
uation: for example m"" for a supercooled liquid will 
be quite different from that for a liquid at higher tem- 
peratures. Third, instead of choosing the scalar current 
densities one could also use each Cartesian component 



recently for a = T, but not for a = i? in order to discuss 
the role of transversal currents on light scattering spectra 
|68| . The resulting equations again are exact but involve 



memory kernels m^'"^/ (<?,*)■ 



B. Linearized equations of motion 

Since the memory kernels will not be independent func- 
tions of the tensorial density correlations S(q, t), the 
third term in Eq. (10b) is a kind of nonlinearity. Within 
MCT m"" (q, t) is approximated by superpositions of 
products S (qi , t) S (q2, t) with q = qi+q2, which makes 
obvious the nonlinear character of the equations of mo- 
tion. It is this nonlinearity which leads to a slowing down 
of the structural relaxation by decreasing the tempera- 
ture or increasing the density of a liquid. This behav- 
ior takes place on the liquid side as a two-step relax- 
ation process characterized by two diverging time scales 
t„ ^\T- Tc\-^ and t - (T - T^y^ where a and 7 are 
positive and Tc is the ideal glass transition temperature. 
The time scale for t„ and r is determined by a micro- 
scope scale to- ^0 depends on inertia and damping effects 
where the latter are due to the regular part of m"" (q, t) , 
accounting for the fast motions. Because our main inter- 
est is the microscopic dynamics, we can get an estimate 
of the microscopic time scale by neglecting the memory 
term including its regular part. This results in a set of lin- 
ear equations from which one immediately obtains for the 
normalized correlator ^{q,t) = 2 (q)S(q, t)S^ 2 (q): 



^{q,t) + n^{q)^{q,t)^0 
with initial conditions 

*(g,0) = l, *((f,0) = O. 
and the hermitian frequency matrix squared: 
n\q) = S-Hq)Y.{q'^frS-Hq) 

a 

or with Eq. (|), (Q and (|l|) 

l"m." 



Irn^Vm" 



M " / 
X (S-^(g) 



m ,i m 



(17) 
(18) 
(19) 

(20) 

(21) 
(22) 



Here, some comments are in order. First, the static 
correlators Sim,i'm'{q) and therefore ^im,i'm'{q) are not 
diagonal in I and V in general. Accordingly, translational 
and rotational modes generally are coupled to each other 
for given q. Second, this coupling vanishes in the limit 
q ^ 0, because the static correlators Sim,Vm'{q) and 
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therefore ^im,i'm'iq) (cf. Eq. ([20|)) become diagonal and 
independent from each other for an isotropic hquid: 



with 



MSo 
IS, 



+ I > 



(23) 

(24a) 
(24b) 



the translational and rotational frequencies for q — > 
and Si = 5;o,io('f = 0). ujo{q) describes the well known 
acoustic isothermal sound wave dispersion and uji (q) the 
optic rotational frequencies for I > 0. 

Third, the limit q was already discussed for a 
molecular liquid using a site-site description |69|. These 
authors also have set the corresponding memory matrix 
to zero, which, by the way, is completely equivalent to a 
short time expansion of the equation of motion in lead- 
ing order. But there the coupling between the partial 
dynamical structure factors does not vanish for q 0. 
In this respect the choice of a molecular representation 
which takes care of the isotropy of the microscopic Hamil- 
tonian is the most natural one, at least for q — > and 
t 0. 

So far we have considered the liquid phase only. The 
idealized version of MCT, which neglects so-called hop- 
ping processes, is a theory developed on the liquid side 
of the glass transition, because it uses as an input static 
correlators in equilibrium. It may also be used below, 
but close to the transition point where the singular be- 
havior still dominates. Whether MCT is even capable to 
describe dynamical feature far below the glass transition 
point is unclear. Nevertheless, we will also apply MCT 
for parameters deeper in the glass phase. An ideal glass is 
a non-ergodic phase with non-ergodicity parameters (not 
normalized) : 



Fem,i'm'{q) — hm S£m,e'm'{q,t) 



(25) 



which are nonzero. One can easily prove that Eqs. (10) 
then imply that 

CZ'rn^' (q) = lini <M (26) 
are nonzero, too, and that 

q"[(C^'''(g1)-T" q"' = (S(91 - F{q))F-\q)S{q)- 

a, a' 

(27) 

Therefore, we introduce S{q,t) and rh"" {q,t) such that 



S(g,t) = F(g) + S(q,t) 
m""'(g,t) = C""' (g) +m""'(g,t) 



(28) 
(29) 



in analogy to the approach in 
Eq. (H) into Eq. (10) yields: 



Substitution of 



S(g,^) -^i^q"S"((f,t) = (30a) 

a. 

S"{q,t)+tq"3''iq)S-\cf)S{q,t) + 
t 

J"(g) J dt'^m""'(<f,t-t')S"'(g,0 + 



+zq"J"(^)S-l(g)F(q) + 
t 

3"iq)J2C°'°''i^ J dt'S°''{q,t') = (30b) 



with initial conditions: 

S(g,0) = S(g)-F(q) 

and 



(31) 



lim S(g, t) = 0, lim m"" (g, t) = 0. (32) 

t — ^oo t — >oo 

With the same argumentation as above we linearize 
Eqs. (^0|b) by taking m"" i'ht) to be zero. Then, taking 
the time derivative of Eq. (pO|b) and using (|30|a) we get 
for S"(g,0 = (J"('?1)-^S"(q,t): 



r(g,i) -f ^(f2'(g))""'S"'(g,i) - 



(33) 



with the frequency matrix squared: 

n\q)-'^' = (J"(g))i [cfS-' (q) cf' - C""'(g)] (j"'(g)) * 

(34) 

Integration of (|30|a) from to i and substituting 
t 

/ dt'S"' {ci^t') from Eqs. (|30|b) with m""' = allows to 



express S{q,t) by S"(g, i): 
S(<f,t) = tF{q)S-' (q) J2 [(C'^"' 



X [r {q)f's^'{q,t) 



(35) 



Similar to fJ, the matrix yCl j is non-diagonal in 

I, I' and in a, a' . This again leads to a coupling between 
the translational and orientational modes, which are the 



eigenmodes of ( Jl ) . Now it is (see M ) 



C"" (q) = q"C"" {q) q" 



(36) 



5 



with: 



limC; 



(37) 



and 



ImA'rn' 'ilnil' 



qT5u'5,nm' (38) 
q, {a,l)^(T,Q) 



1, (a,0^(T,0), 



(39) 



not to be confused with q". Therefore, we get from 
Eq. (H) with Eqs. (i, (0), (H) and (H): 



/ 00,00 



(S"' (^Doo.oo 



M 

l+(s(g)C^^(g) 

^2, 



(40) 



00.00 



for the acoustic part of )"" j. Taking into account 

that we have not normahzed C"" and C"" , the result 
Eq. ( ^ ) is completely analogous to the result derived in 

I^tI for ^ri {q)j for simple one-component liquids. How- 
ever, we note that in contrast to the equations of 
motion for the rescaled correlators S{q, t) and ih"" (g, t) 
are not covariant, due to the splitting of the current den- 
sity into a translational and rotational part. Therefore 

I (g) ) is not a straight forward generalization of 

- 2 

Q, {q) for simple liquids to molecular liquids. 

The coupling between translational and orientational 
modes in the liquid and in the glass already on the linear 
level of the equations of motion is not surprising due to 
the interaction between translational and orientational 
degrees of freedom. But, it is also obvious that mem- 
ory effects will lead to additional couplings and it is this 
point, which we will investigate in section III. 



C. Mode coupling theory 

In the preceding subsection we have neglected the 
memory kernels. Approaching the glass transition sig- 
nificant memory effects occur. Therefore m"" must be 
taken into account. Using mode-coupling theory (MCT) 
an approximate expression for the slow part m"" (q, t) 
has been derived which leads to a closed set of equations 
for S(g, t). This has been done for molecular systems us- 
ing the molecular representation l^j-Q and a site-site de- 
scription |67t|. For a liquid of linear molecules we will 



'iq,t) 



V lala' 

knT 



,{q)m 



hn^l'rn 



slow 



(41) 
(42) 



where the first and second term on the r.h.s. of Eq. 
( ^ ) accounts for the fast and the slow contributions, 
respectively. MCT yields §] 



''lm.,1' 



'{q,t) 



slow 
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(43) 



X VhT.i'rn' {q\qihmi , l[m[ ; g2?2TO2 , ^2"^2) (44) 

X Si^„iul[m[ {qi , t)Sl2m^^l'^m'^ (^2,0 (45) 

The explicit expressions for the vertices V are given in 
Ref. 1^. Equations (10) together with equations ( pl| ) and 
( p3| ) are a closed set of equations for Sim.i'm' {q, t) which 
need the damping coefficients F"^" i,^^, {q) and the static 
correlators Sim,i'm'{q) as an input. The latter uniquely 
determine the vertices. 

It is obvious that the MCT-polynomial ) 

V ' / slow 

leads to additional coupling between the correlators 
Sim,i'm'{q,t). Due to this nonlinearity the MCT- 
equations (10), (|4^) and ( |4^) can only be solved nu- 
merically. We have performed such a numerical solution 
for dipolar hard spheres, one of the simplest systems in- 
volving translational and rotational motion. 



D. Dipolar hard sphere liquid 

The investigation of dipolar hard spheres (DHS) has 
the advantage that approximate analytical expressions 
fZOt for the static correlators are known. 

We consider a system of N hard spheres with homoge- 
neous number density p, diameter d , mass M, moment 
of inertia / = j^AIcP and dipolar moment fi. The ori- 
gin of the body fixed frame is chosen to coincide with 
the center of mass of each sphere, which is the natural 
choice. Since the density of the particles is homogeneous, 
the center of mass is equal to the center of the spheres. 
The reader should note that the MCT-equations (10), 
( ^ ) and (^) and even the original exact equations of 
motion Eqs. (10) are not covariant (i.e. invariant in 
their form) under a shift of the reference point for the 
body-fixed frame. In order to get covariant equations 
one has to project on the individual Cartesian compo- 
nents of J^{q) and j^{q)- Nevertheless, we think that 
the present equations are a reasonable good approach to 
the dynamics of molecular liquids, because of the natu- 
ral choice of the reference frame. Of course, this point 
needs additional investigations. The advantage of MCT 
in site-site representation is that this problem does 
not occur, because no reference point must be chosen. 

The physical control parameters are the packing frac- 
tion (f) — ^pd"^ and the temperature T. In the following 
the length unit is chosen such that d = 1. In addition we 
choose M = 1 and fj, = 1. This choice means that time t 
and temperature T are measured in units of M-^/^d^^'^ / fj, 
and kg^p"^ /d^, respectively. In the following we will use 
T* = T / {p? /ksd^) as dimensionless temperature. 
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As already stated, the MCT-equations (10), (41) and 
(H) require Tf^i,^,{q) and Sira,i'm'{q) as an input. 
Throughout this paper we wiU put all damping coeffi- 
cients r^;,^, to zero. The static correlators Sim,i'm'{q) 
are obtained from Wertheim's solution who used Percus- 
Yevick and mean spherical approximation This leads 
to 



The memory kernels (cf. Eq. (Esh) in the g- frame 



couphngs [||: 



{q,t)6wSmm' contain the following 



Sao{qi,t)Soa{q2,t) and S*!™! (gi, t)5i„i2 (92, t) 

(50) 



S, 



Im.l'' 



'{q) 



Sim{q)Sii'S 

n 

Sll'Smm' 



/ = 0,1, r = 0,1 

otherwise 



(46) 



The vertices of the memory kernels are bilinear in the 
direct correlation functions 



C/m (q) 



47r 
P 



1 - 



1 



Slniiq) 



(47) 



Since we have chosen the so-called q-frame [£T[ (for 
which q — qQ = (0, 0, q)) the correlators and the kernels 
become diagonal in m and m' |^]. The fact that the 
static correlators are diagonal in I and V and are struc- 
tureless for Z, /' > 2 is an artefact of both approxima- 
tions. Additional shortcomings are (i) the independence 
of the center of mass correlator 6*00 {q) on temperature 
and (ii) the smooth behavior of Sim (?) at 0rcp — 0.64, 
the value for random close packing. The static correla- 
tors from Wertheim's approach exhibit a divergency at 
= '/'max = 1- only. We will come back to these points 
below. 

Due to (|6|) and ( ^ only mode coupling occurs be- 
tween the time dependent correlators with Ii,l[,l2,l2 
smaller or equal to one. In order to simplify the MCT- 
equations as much as possible we use the additional ap- 
proximations (in the g-frame) 



Inijl'm' 



iq,t) 



mZpm'iq,t) 



n 





mZ {q,t)Sw5„ 




/ = 0,1, r 

otherwise 



0,1 



(48) 



1^0,1,1' = 0,1 
otherwise 

(49) 



In section IIB we have shown that translational and ori- 
entational dynamics is already coupled on the linear level 
due to non-diagonality in I and Above we have cho- 
sen all correlators to be diagonal in I and I'. Therefore, 
a coupling between Soo{q,t) and Simi,q,t) with I > 
can only originate from mode-coupling effects. There- 
fore, our diagonalization approximation (which is even 
exact in our case due to the restriction of I and /' to 
and 1) allows to study the influence of the mode-coupling 
terms on the microscopic dynamics without interfering 
with the direct coupling mechanism between the correla- 
tors Sim{q,t) discussed in section IIB. 



Wim (<?,*) ^ Sooiqi,t)Sim2{q2,t) and 5i„ii (gi, i)S'oo(<72, t) 

(51) 

Eq. ( pO| ) shows that the center of mass correlator Soo{q, t) 
may undergo a glass transition independently from the 
dipoles whereas Eq. ( [sl] ) demonstrates that the dipoles 
are "slaved" by the center of mass dynamics and can 
freeze only if 500(9, i) has become non-ergodic. In or- 
der that SQo{q,t) and 5i„i(g,i) freeze simultaneously 
the vertices of the bilinear terms in Eq. ( pT| ) must be 
large enough. This happens at low enough temperatures. 
Since Si^-m{q,t) = Si^m{q,t) there are three inde- 
pendent correlators which are the center of mass corre- 
lation function Soo{q,t) and two dipolar ones Simiq,t), 
m = 0,1. In the following we consider the normalized 
correlators (j)im{q,t) = Sira{q,t) / Shn{q)- 

The three static correlators are shown in Fig. |^ for 
two different pairs of [(j),T*). Three main features can 
be seen. First, the g-variation of 6*00 of the "longitu- 
dinal " dipolar correlator and 5*10 is rather similar to 
each other and resembles that of a simple liquid with 
a well pronounced main peak at gmax ~ 27r/a, where 
a Ki d = 1, the mean distance between nearest neigh- 
bors. The "transversal" dipolar correlator Sn [q) be- 
haves quite differently. It exhibits only one peak which is 
at g = and it becomes almost structureless for g > Ag, 
the peak width at half maximum. Second, (g = 0) 
increases with decreasing temperature, whereas Ag de- 
creases. This behavior signals a medium range orienta- 
tional order due to a precursor of ferrofluid order. This 
order is induced by an enhancement of the dipolar inter- 
actions at lower temperatures with respect to the hard 
core repulsion. Third, Sim (g = 0) depends on m, be- 
cause of the long range nature of the dipolar interactions. 

Besides the truncation at I = 1 we must also truncate 
and discretize the g -variable. Since the reduction (for 
— > 00) to a single integral of the sum over gi and 
g2 in Eq. (|3|) for simple liquids ||7|] cannot easily be 
used for molecular liquids the number of steps to calcu- 
late this sum increases quadratically with the number of 
g-values, instead of a linear increase for simple liquids. 
This fact makes the numerical solution of the molecular 
MCT-equations rather CPU-time consuming. Therefore 
we decided to choose a non-equidistant distribution of 30 
g-values between qiow — 0.51 and g„p ~ 40. These values 
were generated by the non-linear relation 



1 



qi, = — arctanh(i^A2,) 



1,2, ...,30. 



(52) 
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with Aa; 



tanh(agc. 



g-^ and Qco — 50. The parameter a has 
been chosen such that the main peak of, e.g. 5*00(9), is 
stiU in the hnear regime of arctanh. The right of Fi g. n| 
demonstrates this for a = 0.065. It is obvious from (|5j) 
that for q < the qi, are almost equidistant and for 
q > a^^ they become more and more diluted. Since the 
large g-regime is not as important as the range around 
the main peak in 5*00(9) or that at q = in 5ii(g) (which 
drive the glass transition) our choice for q^, should not in- 
fluence our results, at least qualitatively. We have solved 
the molecular MCT-equations ( 10), (^) and (||) in time 
space using an algorithm already developed to solve the 
MCT-equations for simple liquids ]73[ |. However, our nu- 
merical procedure differs from that by Gotze et al. for 
the single dumbbell in an isotropic liquid. These authors 
introduced an effective memory kernel (independent on 
a and a') and an effective microscopic frequency. 

Before we come to the dynamics let us shortly discuss 
the phase diagram for a glass transition which was al- 
ready calculated . Throughout the rest of this paper 
all the results are given in the g-frame. The modification 
with respect to Refs. |^,|7^ of the g-discretization leads to 
a small quantitative change of the glass transition lines, 
but without changing the topology. Therefore we have 
calculated again non-ergodicity parameters 



flrn{q) = lim (t)lrn{q,t) 
t — >oc 



(53) 



as a function of (j) and T* from which the phase diagram is 



obtained (Figure El) . Since (f> can not exceed > 



0.64, 



the value for random close packing of hard spheres, we 
have plotted only < (prep- There are two signifi- 
cant temperatures Tj* and Tj. At Tj* the critical line 
'/'type A (^*) (dashed line in Figure ||) reaches (j) = 4>rpc 
and at it merges into the critical fine 0typo bC^*) (solid 
line in Figure ^) . These two temperatures have the fol- 
lowing meaning. For T* > T* an increase of (j) leads at 
•^typo B C^*) to a glass transition for the center of mass 
motion, but not for the dipoles. Choosing T* between 
T* and Tg* again a glass transition of the center of mass 
motion occurs at ^typo b(-^*)- increasing (j) beyond 
(f>type B (^* ) ^ spin glass like transition for the dipoles 
^(T*). Below T2*, center of mass and 

(T*). 



occurs at 



"typo 



%pc B 



dipolar dynamics freeze simultaneously at 
For T* > T2 the critical value (p^y^^ b(^*) is identical to 
that for hard spheres 0^5. For the present choice of q- 
discretization we find (j^HS — 0-5265 which does not differ 
much from the more precise value — 0.516 obtained 
with 100 equidistant g- values |76|| . 

The glass transition at high and low temperatures is 
driven by different physical mechanisms. At high tem- 
peratures it is the so-called cage effect, which leads to 
the freezing of the liquid into a non-ergodic phase, due 
to an increase with increasing density of the main peak 
of the center of mass correlator 5oo (q)- Lowering the 
temperature enhances the role of the dipolar interactions 
which leads to a strong increase of the peak at g = of 



the "transversal" dipolar correlator Su {q). The increase 
of the dipolar correlations are accompanied by a decrease 
of the center of mass correlations (cf. Fig. |l|), and take 
over the role of the cage effect. This behavior refiects 
itself in the non-ergodicity parameters, /n {q = 0) in- 
creases with decreasing temperature much stronger than 
/oo (9max)- For details the reader is referred to [||. 

The properties of the phase diagram require some more 
comments: (i) the T*-independence of •/'typo b(^*) i^or 
T* > Ti originates from the T*-independence of 5oo (q), 
which is an artefact of Wertheim's solution. Removing 
this artefact will result in a r*-dependence of </>typc b (^) 
for all T* with Utiit'^oc ^typo b(^*) ~ '^hS' since dipo- 
lar interactions will be irrelevant at infinite temperature 
(ii), the existence of Tj* relies on the fact that (/'type a(^*) 
approaches the value (/)rcp at a, finite temperature. This is 
true when the static correlators within Percus-Yevick ap- 
proximation are used. Whether or not an improved the- 
ory reproducing the singular behavior at (jj^cp would lead 
to the same conclusion is not obvious. Since such a the- 
ory does not exist, we can not exclude that 0typo a(^*) 
only approaches (prep at T* = 00. 

Now we turn to the time or frequency dependent fea- 
tures of dipolar hard spheres. We have studied the con- 
trol parameter dependence of the dynamics along the 
paths A, B, C and D, indicated in Figure |^. In contrast to 
the determination of the long time behavior inertia effects 
will play an important role in the microscopic time or fre- 
quency regime. In the liquid, these inertia effects enter 
through the microscopic frequency matrix fl{q}. Due to 
the diagonality of S{q) the matrix fl{q) becomes diagonal 
with diagonal elements (in g-frame) which follows from 

Eq. My. 



l{l + l] 



/Sun{q) , (54) 



restricted io I = 0,1. Figure ^ depicts flim{q) for 
two different pairs {4>,T*) in vicinity of the critical line 
"/■typo B (^* ) • This Figure reveals two features which stem 
from the properties of Simil) discussed above. First, 
^im{q = 0) depends on m. Second, and most important, 
a "softening" of the "optic" frequency ilii(g) with q ^ 
is observed for decreasing temperature T* for which the 
role of dipolar interactions (compared to the hard core 
repulsion) becomes more and more enhanced. This "soft- 
ening" comes from the strong increase of 5ii(g) for g — > 
(see Fig. |^) and has its physical origin in the occurrence 
of a medium range orientational order of the dipoles. The 
inverse width (Ag)~^ of the peak in 5ii(g) at g = is a 
measure for the length scale for this orientational order. 
(Ag)~^ increases with decreasing temperature (see Fig. 
nh. We will come back to this point in the final section. 

The microscopic frequencies in the glass phase 
are given by the "renormalized" frequency matrix 

Jl (g) ) I . Taking again into account the diago- 
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nalization with respect to / and I' we obtain from Eq. (^) 
with Eqs. (0) and (||) 



1 



qTq?' +qTqfCZ' (q) Sim (q) 



(55) 



Since C^q" 



(q) is different from zero for a = 
only, we get with Eqs. (i), (pTl), <M) and (|| ): 



noo (q) = f^oo (q) I (1 - /oo [q)) 
C^m iq) is nonzero for aU (a, a') 



1/2 



a' = T, 



(56) 



Therefore two 

eigenfrequencies f2]^„j (g) exist for each m. tloo {q) and 
^tm il) shown in Fig. |4| the "renormahzed" fre- 
quencies require the non-ergodicity parameters as an in- 
put. The latter become less accurate for small q, be- 
cause the discretization of the q-values influences the 
results much stronger for decreasing q. Therefore we 
do not present data below q = 1. Of course it is 
^oo{q) cq and (if^iq) ^ tlf^iO) > for q ^ 0. 
Note that (i) Cl^-^(q) does not vary much with q and (ii) 

that ^iiiq) < ^ooiqmax ^ 7). 

Now let us turn to the equations (10) including the 
time-dependent memory kernels given by ( ^l|) and (43). 
From the solutions of these MCT-equations we get the 
normalized correlators (j)im {q, t) or equivalently the cor- 
responding susceptibility spectra x'lmil^ related to the 
correlation spectra (/)"„(g,[j) by xL(9,w) = uj(l>[[^{q,uj). 
Figures || and | show </>;,„ (q, t) and x'/^iq, to) for q ~ 4.7 
along path A in Figure ||. Figure || clearly demonstrates 
for all correlators the formation of a plateau and the 
tremendous slowing down of the relaxation with increas- 
ing which stops at (j)c{T*) ~ 0.3782. 

The dynamical features are better recognized in uj- 
space. Before we discuss these results, let us comment on 
the quality of our w-dependent data. Because the glassy 
dynamics involves many decades in time one has to use 
a special algorithm, which decimates the time step with 
increasing time, i.e. one chooses an initial time increment 
h to discretize the time axis [ |73| . After a certain number 
of time steps the increment is doubled, etc. Although 
the i— dependent data shown in Fig. ^ look continuous, 
they exhibit tiny discontinuities of the correlators and 
of their first time derivative at those times at which the 
increment is doubled. These discontinuities which could 
be diminished by an increase of the CPU-time lead to 
the wiggles which can, e.g. be seen in Fig. ^ However, 
we have checked by variation of the parameters h, the 
boundary values of the Fourier integrals, etc. that the 
features we will address below are not influenced by the 
wiggles. Further, the low frequency wing of the a— peak 
for x'lmi'i^^) I^ig- H shows a rather small deviation 
from a linear w-dependence, as it should be. This failure 
is due to the choice of the lower bound for the Fourier 



integral. If this is taken to be zero, the wiggles become 
more pronounced. Therefore, choosing a nonzero bound 
is a compromise between the reduction of the wiggles and 
a small deviation from the linearity in uj for w — > 0. We 
have also used a spline-technique for smoothing the data 
in time space. The Fourier transform of those has repro- 
duced all the relevant features obtained from the original 
i— dependent correlators. 

Let us now discuss the spectra. For (f> < (j)c{T\), i.e. 
in the liquid phase, Figure left column reveals the exis- 
tence of three peaks for Xooiq — 4.7, uj) and two peaks for 
x'liiq — 4.7, w) and Xioil — 4.7, cj). The low frequency 
peak in all spectra x'l'mil — 4-7, w) is the a-peak related 
to the slowing down of the structural relaxation (transla- 
tional and rotational degrees of freedom). Xooil — 4-7, w) 
has a high-frequency (hf) peak at ujhf ~ 10 and an ad- 
ditional peak at ujop — 1, i.e. about one decade below 
the hf-peak. x'ioiq — 4.7, w) and Xii('Z — 4.7, w) ex- 
hibits only one peak in the microscopic frequency do- 
main at cj « LOop- Below we will show that the peak in 
Xoaiq — 4.7,0;) at w « Wop originates from the orien- 
tational dynamics. Therefore we call it " orientational" 
peak (op). 

The situation for (p > (f)c{T^), i.e. in the glass phase, 
is similar with two exceptions (cf. right column of Fig. 

First, it is clear that there is no a-peak anymore and 
second an additional peak located between the op and 
the hf-peak exists in Xoa{q — 4.7, w) for (f) ~ 0.53 and 

0. 60, i.e. deep in the glass phase. Even a further "peak" 
looking more like a shoulder appears for (j) = 0.6. Such 
additional peaks can also be seen in xi'oCs — 4-7, w) and 
Xiiiq — 4.7, w) (cf. right column of Fig. ^) 

In the following we will concentrate ourselves onto 
the microscopic regime. Therefore we can use a linear 
tj— scale, which allows us to recognize the u- and (f>- 
dependence (which we have described above) much bet- 
ter for 0.1 < < 10. This is done in Figure for the 
data from Figure ^ Both microscopic peaks at uJop ~ 1 
and ujfif ~ 10 can clearly be seen in the compressibil- 
ity spectrum Xooiq — 4.7, w). Whereas the intensity of 
the hf-peak increases with increasing 0, the op becomes 
more pronounced when approaching the glass transition 
from the liquid side. In the glass phase at </> = 0.44 and 
4> = 0.48 it is less prominent but again becomes more 
pronounced for (j) > 0.5, although its intensity decreases. 
The appearance of the intermediate peak for (j) = 0.53 
and (f) — 0.60 at Lo K 6 and lo ^ 8, respectively, can be 
observed as well as the shoulder at a; « 4. Due to the 
linear w-scale we clearly see that Xioil — 4-7, o^) also ex- 
hibits an intermediate peak for </> = 0.53 and = 0.60, 

1. e. Xioil — 4-7, w) resembles Xaoiq — 4.7, w) but with 
the opposite ^-dependence of the intensity of the op and 
hf-peak for (j) > (j)c{T^). x'iiiq — 4.7, w) which is propor- 
tional to the dielectric loss e"{q ~ 4.7, w), possesses only 
one well pronounced peak at lo k ujop- To check how far 
the features depend on the path through a critical point, 
we have investigated x'lm^q — 4.7, w) along path C (see 
Figure |l|) for which (p = 0.3786 is fixed. This allows 
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to study how far density and temperature variation lead 
to similar or different susceptibility spectra. The results 
for path C are given in Figure||. x'ooil — 4-7, uj) shows an 
op and a hf-peak at tOop ~ 1 and LOhf ~ 10, respectively. 
Their corresponding intensities behave qualitatively sim- 
ilar for decreasing T* as we have found along path A for 
increasing 0. Intermediate peaks between op and hf-peak 
do not occur. x'loiQ — 4-7, uj) differs from the result along 
path A, whereas Xii{q — 4.7, looks similar. 

Translational and orientational dynamics freeze on 
path A and C simultaneously. Therefore we have also 
calculated x'lmi^ — 4.7, a;) on path B for fixed T* = 
0.30 where 4>m{l,i) freezes at cf)^^ ~ 0.5265 first, 
whereas (l)im{<],t) undergoes a spin-glass-like transition 
at ^type a(^* = 0-3) — 0.62. The corresponding results 
are presented in Figure ^ Xooiq — 4.7, w) and x'lmil — 
4.7,0;) exhibit a well pronounced main peak. Its position 
depends sensitively on in case of x'mil — 4-7, and is 
practically (j) independent for Ximi.1 — 4.7,0;). The lat- 
ter quantity also pocesses a small peak at uj w 10, which 
originates from the translational motion via translation- 
rotation-coupling. 

A peak in the susceptibility spectra can be either of 
relaxational or oscillatory origin. To distinguish between 
both types of behavior one must study the o;-dependence 
of the correlation spectra 0"^(q,o;). Spectra from neu- 
tron scattering are superpositions of 0"^(g,o;) |^^. A 
peak in (j)'l^(q,ijj) at o; > proves the existence of an 
oscillation whereas a peak at o; = is of purely relax- 
ational type. The corresponding peak width is a measure 
of the damping. These correlation spectra are shown in 
Figure |l^ and |ll| for path A and path B, respectively. 
For (j) < (t>ciT*) (not shown in Fig. |l^ and 11) there 
is no peak at nonzero frequency, except for the hf-peak. 
Therefore the op we have found in Ximil — 4-7, o;) in the 
liquid phase is a purely relaxational excitation. Now let 
us discuss the correlation spectra in the glass phase. Fig- 
ure |l^ shows that in the glass, but close to the critical 
packing fraction, there is a hf-peak at o)/i/ ~ 10. Deeper 
in the glass the position of that hf-peak shifts to higher 
frequencies and an op at u>op ~ 1 appears. Decreasing 
even more an intermediate peak between the op and the 
hf-peak is produced as well. The position of the op and 
intermediate peak shifts to higher frequencies with de- 
creasing similar to the hf-peak, due to an increase of the 
glass-stiffness. The 0-dependence of 4>'imil — 4.7, o;) at 
much higher temperature T* = 0.3 (path B in Figure |^) 
is presented in Figure |l^. We observe a similar behavior 
as in Figure |l^, i.e. besides the hf-peak at tOhf = 20 — 60 
for 0.56 < 4) < 0.63 there occurs the op at Uop ~ 1 and 
for (f) = 0.60 and = 0.63 a shoulder at about o; « 20 
which corresponds to the intermediate peak in Figure 

One of the conclusions we can draw from x'lmil — 
4.7,0;) and 4'imil — 4-7, o;) is that there is an op at 
LUop ~ uJop ~ 1 roughly one decade below a hf-peak. The 
hf-peak relates to a damped oscillation for the investi- 
gated range of T* and (j), whereas the op changes its char- 



acter from purely relaxational to damped oscillational 
behavior under going from the liquid to the glass. Since 
the dipolar spectra Ximil — 4-7, o;) and (f>ijn{q — 4.7, o;) 
have a main peak at about o; « 1 it is tempting to asso- 
ciate the op in Xooil — 4.7, o;) and 000(9 — 4.7, o;) with 
the orientational degrees of freedom and their coupling 
(via mode coupling effects) to the translational ones. If 
this interpretation is correct then the op must exhibit 
an isotope effect with respect to a change of the mo- 
ment of inertia /. Figure |2| presents x'lmil — 4-7, o;) at 
TX = 0.04, (j) = 0.381 for / = (the correct value for a 
hard sphere with M=l and d=l and homogeneous mass 
distribution), 1 and 10. From this result we find 



1 



OJop OC 



(57) 



This type of isotope effect also occurs for 4>'i'„^{q, o;) which 
is shown in Figure |l| for (/);'„^(q ~ 4.7, o;) at = 0.04, 
(/) = 0.530 and for / = j^,! and 10. This result yields: 



(58) 



The approximate scaling of Wop and iliop with 1/VT 
strongly supports the orientational origin of the op. Since 
LUfi f and LUhf are rather insensitive on a change of / their 
origin must lie in the translational motion. 

So far we have not studied the g-dependence of all of 
these peaks. Since this may elucidate more the features 
of the several peaks, we present x"^(g,o;) at T4 = 0.04, 
(f) = 0.381 in the glass but near the glass transition and 
(/)"„((7,o;)at TX = 0.04, = 0.53 deeper in the glass in 
Figures and respectively. The result for Xooi^^^) 
shows a nearly (/-independent position of the op at oj^p « 
1 for l.O < q < 10.6, whereas the position oj^j of the 
hf-peak changes with q. A similar q -sensitivity holds for 
the position of the main microscopic peak in Ximily ^) fo'^ 
TO = and 1. The result for (j)'i^(q,uj) yields the same q- 
independence of the op at ojop — 1 and a high sensitivity 
of the position of the hf-peak at dihf ~ 10. The high 
frequency peak is practically absent within a "window" 
around q = 6.5 and it appears below and above that 
window. 

Finally, let us comment on the dependence of the peaks 
on T* and (j) as well as on a, which characterizes the dis- 
tribution of the values for q (cf. Eq. ( [5^ ) and Figure]^). 
Let us start with the a-dependence. Figure |l^ presents 
XooiQmax.uj) at = 0.04, = 0.381 for three different 
values of a. qmax which depends on a has been chosen 
as that g- value closest to the main maximum of S'oo('?). 
Although the peak positions and also the intensity varies 
with a the qualitative features do not depend on a, at 
least for reasonable chosen values for a. 

The T*-dependence of the position ujop and of the 
height hop of the op for q ~ 4.7 is shown in Figure 
|l7| a for 4>£i — 0.525 along path D. The corresponding 
(f) dependence for the same q value and TX = 0.04 along 
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path A is given in Figure p7|b. hop follows a linear T* 



fa = lim 0Q(t). 



(63) 



dependence between T* 
we have studied) and T* 



0.04 (the lowest temperature 
0.1. In a temperature region 



where ujop has a minimum, hop shows a crossover to a 



constant. As can be observed from Figure 173 the po- 
sition LOop and the height hop are almost constant below 
4'ciTX = 0.04) ~ 0.3782 and increase with above that 
value. 



III. SCHEMATIC MODEL 

In this section we will demonstrate that besides a 
high frequency peak the occurrence of an additional mi- 
croscopic low frequency peak can be obtained from a 
schematic model. For dipolar hard spheres it has been 
shown in section |ll| that this additional peak in Xooil^ 
and (f)QQ{q,uj) has its origin in the orientational degrees 
of freedom and their coupling to the center of mass mo- 
tion. Therefore the schematic model must contain at 
least two correlators 4>Q{t) and <j)i{t), describing center 
of mass and small q orientational dynamics, respectively. 
The topology of the phase diagram (Fig. ||) is related 
to the special form of the memory kernels (Eq. ( |50| ) and 
(|l|)), i.e. TO^o"' does not only involve Soo{qi,t)Sooiq2,t) 
but contains also a feedback from the dipolar motion, 
due to Simi iqi,t)Sim2 (92, t). In contrast m"^ {q, t) only 
possesses a coupling between the center of mass and the 
dipolar correlators, i.e. the latter can only become non- 
ergodic if the center of mass correlator becomes frozen. 
A schematic model which exactly describes this behavior 
is the Bosse-Krieger model 



(59) 



JQ 



a = 0, 1 with the memory kernels 

mo{t) = Fo(0o(i),0i(O) =a0oW+6</'?(O(6Oa) 
mi(i) = Fi{Mt),Mt)) ^ ^3Mt)Mt) (60b) 

and initial conditions 0a (0) — 1, (j^ait) — 0. To mimic 
the (j)- and T-dependence of the dipolar hard spheres we 
choose 



^1 



6 = x/T, 6 = l/T. 



(61) 



This choice is motivated by the observation that the 
mode coupling between center of mass and dipolar dy- 
namics vanishes for T — > 00, leaving only the cage effect 
for simple hard spheres with "packing fraction" 0. x is 
considered as an additional parameter. From ( |59| ) and 
(60) we immediately get the equations 



fa 



FaifoJi), a =1,2 



1-fa 

for the corresponding non-ergodicity parameters 



(62) 



( |6^ ) has been solved for x = 0.15. From this solution we 
get the phase diagram for glass transitions presented in 
Figure Comparison with Figure ^ reveals qualitative 
agreement between both phase diagrams. We emphasize 
that this result holds only in a certain range for x. The 
full Ci ~ ~ $3 or equivalently (f> — T — x diagram was 
explored in detail in Ref. Q. Since /i = in the glass 
phase II (see Figure Eqs. ( |62| ) for a = reduce to the 
.Fj-model |I] for which the critical value for (p is = 4 
and the critical non-ergodicity parameter is f^' = ^ (see 
horizontal line in Figure p^ . 

Having found the same qualitative phase diagram as 
for dipolar hard spheres we can now investigate the dy- 
namics for the Bosse-Krieger model. This has already 
been done in order to calculate the long time dynamics 
for the a-relaxation and the a- and relaxation . 
An investigation which also includes the microscopic fre- 
quency regime can be found in Ref. Q . Figure 4b of this 
Reference gives a hint for the existence of an extra peak 
between the a- and the hf-peak. But this phenomenon 
has not been discussed there. 

In Figure |l^ we present (f>o{t) for fixed frequency flo = 
1 and different values for f2i. For f2i = f2o we find the 
typical two step relaxation process from (/'o(O) = 1 to 
a plateau at fc = f^'^ and a final decay to zero. With 
decreasing fli we observe that a shoulder occurs which 
finally develops into a second plateau above f^^ This can 
easily be understood. In case that becomes more and 
more soft, with respect to Qq, the second term in mo(t) 
(cf. Eq. (60a) varies slower and slower, already on the 
microscopic time scale to = ■ Therefore the relax- 
ation kernel mo{t) has a rather slowly varying part. It 
is this part which generates the second plateau. Since 
the decay from a pronounced plateau produces a relax- 
ation peak in the corresponding susceptibility spectrum 
we expect, besides the high frequency peak related to Qq 
and the a-peak, one more peak. This peak stems from 
the decay from the upper plateau to that at f^^. Fig- 
ures |o| and ^ present Xa(^) ^0 = Ij ^^i = 0.1 and 
1^0 = 1^1 = along path A' and B' (see Figure respec- 
tively. Starting in the ergodic phase not too close to the 
transition line there is in Xq(w) the a-peak (visible on 
the linear scale for e < and n = 1, only) and the high 
frequency peak at w w fio, only. Approaching the criti- 
cal line, we observe that the high frequency peak shifts 
to higher frequencies, due to "stiffening" (see Ref. |^), 
and the a-peak to lower frequencies. With this shift of 
the a-peak, a third peak becomes more pronounced at a 
position which is about the same as the position of the 
main peak in Xii^) and which scales with the bare fre- 
quency ill. This peak remains also in the glass phase. 
To investigate the nature of that additional peak again 
we have plotted (t>ai^) in Figure |2^ and |2^ for path A' 
and B', respectively. The hf-peak occurs already in the 
"liquid" phase (at least for the lower "temperature") in- 
dicating a damped oscillatory dynamics. The spectral 
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shape of the hf-peak strongly resembles the shape of 
the anomalous oscillation peak (AOP) found in Ref. ||5^] 
for a hard sphere system, a one-and a two component 
schematic model. Therefore our notation hf-peak should 
not be confused with the high-frequency sound peak in 
Rcf. In contrast, the additional peak at a; « fii does 
not exist in the "liquid" but in the "glass" phase, only. 
Hence, the character of the extra peak changes from a 
relaxational to a damped oscillational type. Therefore 
the occurrence and the qualitative features of the addi- 
tional peak are quite similar to the behavior found for the 
system of dipolar hard spheres. Since Qi mimics the ori- 
entational frequency nim{q), we have further support of 
the orientational origin of the orientational peak found 
in XonilT^) ^-nd for dipolar hard spheres. In 

addition, the schematic model also demonstrates that an 
additional microscopic low frequency peak can only exist 
when the microscopic dynamics for 01 (t) is "soft" enough 
with respect to that of 0o(O- 

For completeness we mention that for the lower "tem- 
perature" (path A') there exists in (j>Q{uj) a further peak 
at a; w 0.4. This peak becomes more prominent on a 
double-logarithmic representation (see inset of Figj2^) 
and its position also scales with 17 1, and not with ^Iq. 
Hence, both peaks at a; w 0.1 and uj w 0.4 have their 
origin in the "orientational" motion of 0i(t) and its cou- 
pling to 4>oit) via mode coupling terms. Performing a 
"stiff-glass" -approximation ||5^ for the BK- model (which 
has not yet been done) should allow to elucidate the peak 
structures and the spectral shapes more in detail. Since 
this is not the purpose of the present paper this will be 
discussed elsewhere. 

It is also instructive to neglect memory effects for ipi (t) , 
i.e. we put mi(i) equal to zero. Then, ( |5^ for a = 1 de- 
scribes a damped harmonic oscillator with a feedback on 
(j)o{t) due to the second term in mo{t). For this "reduced" 
Bosse-Krieger model we present x'oi'^) Figure ^ for 
= 1, ill = 0.1, 1^0 = and i/i = 0.1. Note that we in- 
troduced a nonzero damping i^i in order to comply with 
the nonzero width of the main peak in Xii^) for tl^^ 
original BK-model. As expected one observes an addi- 
tional peak which originates from the damped harmonic 
motion of (/'i(t). Such an observation has already been 
made earlier for another type of schematic model . 
Comparing Xoi^) from Figure |2l| and (/'o(w) from Fig- 
ure with Xoi'^) ^-n^d 0o(w) in Figure ^ shows that 
the memory effects for (/'i(i) do not influence Xoi^) s-nd 
4>Q{Lij) qualitatively. 



IV. DISCUSSION AND CONCLUSIONS 

Recently we have extended the mode coupling theory 
(MCT) for simple or binary ||l|-Q to molecular liquids 
P,[7pq| . This has been done within a molecular represen- 
tation, which separates translational from orientational 
degrees of freedom. This molecular theory has already 



been applied and tested for diatomic molecules || 8 ^, p^ , p^ 
and water molecules with SPC/E potential |]7|,|l4|. These 
tests were restricted to the comparison of the non- 
ergodicity parameters and critical amplitudes (for di- 
atomic molecules, only) from molecular MCT with the 
corresponding quantities from a MD-simulation. A sat- 
isfactory agreement between MCT and MD-simulation 
has been found So far no dynamical results were deter- 
mined from mode coupling theory for molecular liquids. 
The only solution of the time-dependent MCT-equations 
for a system with orientational degrees of freedom has 
been obtained recently for the orientational correlators 
of a single dumbell in an isotropic Hquid |^,^. In the 
present paper we have solved the time-dependent molec- 
ular MCT equations for a system of dipolar hard spheres. 
This is one of the simplest systems involving translational 
and orientational degrees of freedom. In addition it has 
the advantage, that the static correlators which are the 
input quantities for the MCT-equations are known an- 
alytically within some approximations. These approx- 
imations make the tensorial correlators in the g-frame 
(t>im,i'm'{q = qo,t) = (t)im{q,t)6ii'6mm' diagonal in I and 
I'. In addition they lead to a restriction of 4)im{q,t) to 
Z = 0, 1 only. This means that our calculations yield 
information on center of mass and dipolar dynamics. 

The major result we have found in the microscopic fre- 
quency regime is the existence of an orientational peak at 
^op (hquid) « cj op (glass) about one decade below the high 
frequency peak. This orientational peak exists above 
and below the glass transition lines in the two dimen- 
sional phase diagram of a dimensionless temperature T* 
and packing fraction provided that the a-peak posi- 
tion oja is much smaller than to op- On the liquid side it 
exists in XQQ(q,w), the compressibility spectrum but not 
in 0QQ(g,ti;) which proves its pure relaxational character. 
However, crossing the glass transition line this peak also 
appears in the correlation spectrum 0QQ(g,u;) if T* and 
(j) becomes, respectively, small and large enough. This 
means that the orientational peak exists on both sides of 
the glass transition lines and changes from a relaxational 
to a damped oscillational type of excitation under trans- 
formation of the liquid to the glass. The manner how the 
op develops in Xoo('?, oj) as a function of 4> (see the double 
log representation in Fig. |^) resembles at least qualita- 
tively, the behavior of the extra peak in the susceptibility 
spectra for e.g. Salol and Orthoterphenyl [ p3| , p0p^ ] at 
Lo K, 300 GHz under variation of the temperature. 

Since the position ujop is almost independent on q it is 
related to a localized, non-propagating mode. In addi- 
tion, the isotope-effect, i.e. tOop oc /~^/^, clearly proves 
its orientational origin. Since the bare frequencies ^i„i (q) 
for Z = 1 scale with /"^/^ for q = 0, only (cf. Eq. (p4|)), 
the orientational peak must be generated by orientational 
modes with q 0. Why do so long- waved orientational 
modes play a crucial role? An answer follows from the 
static correlator Sii{q). In contrast to Sooiq) and Sio{q) 
it exhibits a dominant peak with width Aq at g = and 
decays rapidly to one for q > Aq (see Fig. |l|). With 
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decreasing temperature its peak height increases and Aq 
decreases. The increase of Sii{q ~ 0) has two imphca- 
tions: 

(i) since the glass transition in MCT is driven by the 
growing of the main peaks of the static correlators R, it is 
the increase of Sii{q ~ 0) and the coupling of the orien- 
tational mode pii(q,t) to pQo(q,t), the translational one, 
which enhances the tendency for glass formation. This 
type of glass transition mechanism was already described 
within the framework of MCT for crystalline systems. 
For orientational glasses I 



and strictly periodic lat- 
tices 88 1 which undergo a second order equilibrium 
phase transition the increase of the critical fluctuations 
already lead to an ideal glass transition above the equilib- 
rium transition temperature. The increase of the critical 
fluctuations is accompanied by the softening of the criti- 
cal mode. 

(ii) since nn{q) cx 1/y^S^qj (cf. Eq. (|| )) the 
increase of Sii{q ~ 0) implies a softening of the orienta- 
tional mode for g « (see also Fig. ^. It is this soft- 
ening of the orientational mode with g « and I — I, 
m — 1 which results in a frequency scale separation be- 
tween this orientational mode and the translational one 

= 0, TO = 0) with q « qmax- Therefore two peaks, the 
orientational and the hf-peak, become visible. The con- 
nection of coopiq) with the (1, l)-mode is also supported 
by the fact that uJopiq) « ^fiil) (see Fig. ||). 

The narrowing of the peak in Sii{q) at g = also has 
an interesting implication. Since l/Aq is a measure of 
a correlation length, the narrowing implies a growth of 
an orientational order. Although l/Aq does not diverge, 
it may be significantly larger than a few Angstrom, i.e 
there may exist a medium range orientational order. For 
instance, for (j> ~ 0.381 and T* — 0.04 the correlation 
length is about six diameters (cf. Figure 0). The fact 
that the increase of Sii{q « 0) is accompanied by the 
narrowing of the peak in Sii{q) at q — proves the ex- 
istence of a correlation between the occurrence of the 
orientational peak and medium range orientational or- 
der. This is an interesting observation insofar that such 
a correlation has already been predicted ||8^- |9l| , |39| , ^ , 
although this does not seem to be a universal feature 
p3| . We emphasize that we do not consider this type of 
behavior as an exceptional case which exists for dipolar 
hard spheres, only. For example for a system of hard el- 
lipsoids which may exhibit medium ranged nematic order 
(see Ref. [Q) similar results are expected. 

Besides the op, we have found in Xooil:^) and 
(j)QQ{q, uj) above and below the glass transition a high fre- 



^That the underlying mechanism leading to an ideal glass 
transition can be much more sophisticated, i.e. not only re- 
lated to the increase of e.g. one peak in 5*00(9) has been 
recently demonstrated f or a system of hard spheres with at- 
tractive interactions |9£-102| 



quency peak at tUhfiq)- It exists in the investigated range 
0.51 < q < 40 except in an interval centered around qmax, 
the position of the main peak in Sooil)- ^hf{q) is shown 
in Figure ^ for q < 2. Despite an inaccuracy of about 
10 (due to the g-discretization) its (/-dependence is in 
phase with the g- variation of the renormalized frequency 
f2oo((7). These qualitative features are similar to the type 
of the high frequency peak found for a system of hard 
spheres Q. Since the accuracy ofujhf{q) becomes even 
worse for g ^ it is not possible to detect the expected 
linear g-dependence for q —^ 0. Accordingly, we can not 
yet prove that the hf-peak for DHS corresponds to the 
high frequency sound, as found for the system of hard 
spheres ||5^. In ||5^ it has also been shown that an ad- 
ditional peak, which was called "anomalous oscillation 
peak" (AOP), appears deep in the glass phase. Its origin 
lies in the harmonic motion of the particles in their cages 
prf . Whether or not the hump between the op and the 
hf-peak we have observed in 0oq((7, cj) deep in the glass 
(see Fig. and p^) corresponds to the AOP is not clear. 
Again due to the restriction to 30 values for q and also a 
smaller number of discrete values for t compared to Ref. 
p?! we can not give detailed quantitative information 
on, e.g. the g-dependence of the position and width of 
that hump. On the other hand the high frequency wing 
of the orientational peak may interfere with the AOP, 
thus complicating the analysis of the latter. 

A well-pronounced two-peak structure in the micro- 
scopic frequency regime can also occur for XiailT^) (cf. 
Fig. 1^) deep in the glass, whereas this does not happen 
for Xii(9-w) (cf. Fig. |^) which is directly related to the 
dielectric loss spectrum e"{q,uj). 

Let us come back to the orientational peak. Figure p^a 
shows that its intensity varies linearly with temperature 
below about T* = 0.13 and for (j)D = 0.525, i.e. for T* 
smaller than T*{<j)D = 0.525) ~ 0.13. This fact, as well 
as the rather insensitivity of its position cuop on q and 
its location about one decade below the high frequency 
peak is reminiscent to the boson peak. Additionally the 
increase of uJop with packing fraction is similar to the shift 
of the boson peak to higher frequency by increasing the 
pressure [p5|-p7t. 

To conclude we can say that we have found in the mi- 
croscopic frequency regime of the compressibility spec- 
trum and the corresponding correlation spectrum an ad- 
ditional orientational low frequency peak. It originates 
from a localized and non-propagating orientational mode 
coupled to the longitudinal acoustic sound waves (trans- 
lational modes). Its dynamical origin changes from a 
pure relaxational to a damped oscillational type when 
crossing the glass transition from the liquid side. Since 
the orientational frequency has a gap at g = the orien- 
tational excitation is an optical mode. These character- 
istic features coincide with several explanations for the 
existence of the boson peak ]l8| , ^ , ^8|j6^j6^ which may 
stress the role of orientational degrees of freedom. It is 
also interesting that the intensity of the op in a certain 
temperature regime (cf. Fig. [l7|a) changes linear with 
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T* which is also true for the boson peak. In addition we 
have found an interesting correlation between the orien- 
tational peak phenomenon and a medium range orienta- 
tional order, as already suggested earlier for the boson 
peak [89|- |9l] , |3^ , |9^ , although it has been shown that this 
does not seem to hold for all glass formers j93| . 

For dipolar hard spheres this order is of ferroelectric 
type and for hard ellipsoids related to a precursor of ne- 
matic order. For Si02 it can occur due to strong orienta- 
tional bond interactions. As demonstrated in the present 
paper this medium range orientational order may result 
in a softening of a localized, optic orientational mode 
and may finally generate an additional microscopic peak 
about one decade below the high frequency peak. 

In general it is not so easy to separate from experimen- 
tal spectra obtained from light or neutron scattering the 
contribution from the orientational motion. But this is 
not true for a MD-simulation. Due to the availability of 
all the microscopic information it would be desirable to 
explore in more detail by MD-simulations the role of the 
orientational degrees of freedom for the glass transition 
itself and also for the spectra in the microscopic regime, 
as it has already been done for the sound propagation in 
liquid water |Q. 

Finally we mention that the microscopic spectral fea- 
tures we have found for DHS can be qualitatively repro- 
duced by use of a schematic model, the Bosse-Krieger 
model. 
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APPENDIX A: THE NUMERICALLY SOLVED 
MCT EQUATIONS 

Using Eq. (^ we get the first equation of motion for 
the generahzed density-density correlation functions (in 
the g-frame). 



Sll'm{q, t) 



a=R,T 



with 



(Al) 



(A2) 



the density - current- density correlation function. Pro- 
jecting on the tensorial densities (Eq. (^) and the lon- 
gitudinal translational and rotational current densities 
(Eq. (5)) at once leads to 



E 



(A5) 



where '^ii^^iq^t) is determined only up to an integra- 
tion constant which does not influence the result for 
Su'm{q,t). Making use of (A4) and taking the time 



derivative of Eq. (AS) leads to a set of equations of sim- 
ilar structure as the MCT equations for simple liquids 
[|l|,||j|J||. The molecular MCT - equations do not only 
couple all (discretized) wave vectors, as in simple liquids, 
but also the indices a — T, R for translation and rota- 
tion and the spherical indices I < Ico and —l<m<l. 
The functions are auxiliary functions because their in- 
troduction is just for numerical purpose. The integral in 
Eq. (A3) transforms to 



dt'MrUq,t--t')'l> 



(A6) 



Now a time derivative is multiplied with the integration 
measure dt. As a conseque nce the time step Ai can- 
cels at a discretized form of (A6). Using the decimation 
technique [^ At becomes very large. In this case the 
calculation of the integrals would become unstable if At 
would not cancel. 

Since the auxiliary functions only occur as time deriva- 
tives in the original equations the result is independent 
on the initial values of ^uij^iq^t = 0). When solving the 
molecular MCT equations we have chosen: 



for I = Oyi' and 







(A7a) 
(A7b) 

(A7c) 
(A7d) 



for I > as initial values. 



knT. 



E qr{q)^is-\q))ihmSi,i'rn{q,t) (A3) 

l2—m 

-EE ^-Y^ f dt'M^,::M.t-t')<pCL^^^^^ 

'■'=ma'=R,T " "^^ 



with the cut-off value Ico for / (for DHS it is Ico — ^, cf. 
Eqs. (H), (H)) and Ic defined in (^. To make use of 
the established numerical method to solve such equations 
[p^ it is necessary to introduce auxiliary functions <i>: 



(A4) 



As a consequence Eq. (Al) reduces to 
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FIG. 1. The q dependence of the static structure factors 
Sim{q) for (Im) — (0,0), (1,0), (1,1) and two different pairs 
of (0,r*) indicated by the full circles in Figure |: (0.53,0.3) 
(left), (0.381,0.04) (right). Aq is the peak width of 5'ii(<3) at 
half maximum. The symbols indicate the discretized g-values 
(see text) 
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FIG. 2. Glass transition phase diagram for dipolar hard 
spheres: The solid and the dashed line are critical lines at 
which a discontinuous (type B) and a continuous (type A) 
glass transition takes place between the phases I, II and III. 
A, B, C and D denote the various paths on which we have in- 
vestigated the control parameter dependence of the dynamics. 
The full circles indicate the two points at which the static cor- 
relators in Figure |l] and the microscopic frequencies in Figure 
^ were calculated. For Ti and T2 see text 
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FIG. 3. g-dependence of the microscopic frequencies 
n,™(q) for {l,m) = (0,0) (solid), (1,0) (dotted) and (1,1) 
(dashed) (a) cj> = 0.381, T* = 0.04 and (b) (j) = 0.53, T* = 0.3, 
indicated in Figure ^ by the full circles 
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FIG. 4. g-dependence of the renormalized microscopic fre- 
quencies for (f) — 0.53, r* = 0.04: (a) f2oo(g) is shown as solid 
line; the position Uhf and ujop obtained from <^oo(g, <^) are pre- 
sented by the full circles and spheres, respectively (b) ^^^{q) 
and fjfj (g) are shown as circles and spheres, respectively. 
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FIG. 5. lin-log plot of (l>im{q,t) at q ~ 4.7, TX = 0.04 and 
(/!> along path A in Figure ^ The dashed line corresponds to 
the results at 6 ~ cbc — 0.3782 
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FIG. 6. log-log plot of xL(g,'^) at g ~ 4.7, TX = 0.04 and 
(j> along path A in Figure E| liquid (left column) and glass 
(right column). The dashed line corresponds to the results at 
(l> — 4>c — 0.3782. (^-values for Xio(9i'^) and Xii{q,^) are the 
same as for Xooilj^) 
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FIG. 7. x"m(?,^) on linear scales at g ~ 4.7, TX = 0.04 
and (j) along path A in Figure ^ The dashed line corresponds 
to the results at <j!> ~ <j!>c — 0.3782 
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FIG. 8. x"m('7)^) on linear scales at g ~ 4.7, (^c = 0.3786 
and T* along path C. The dashed line corresponds to the 
result at T* ~ T* ~ 0.04 
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CO 

FIG. 9. x"m(9>^) on linear scales at q ~ 4.7, Tg = 0.3 
and </> along path B. The dashed line corresponds to the result 
at ~ <^c — 0.5265. Because Ximil,^) does not vary much 
with we have not labelled the various curves with (j) 




FIG. 10. (p'l'miQ — 4.7, oj) on linear scales in the glass phase 
along path A. I.e. T* = 0.04 and <^ = 0.381 (solid), = 0.44 
(dotted), (j> = 0.48 (dashed), (j) = 0.53 (long-dashed) and 
4> = 0.6 (dot-dashed) 
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FIG. 12. x'Lig - 4.7, tj) at <;/> = 0.381, T 
FIG. 11. (f>'L{q ~ 4.7, w) on linear scales for T* = 0.30 and / = _^ (solid), 1 (dotted) and 10 (dashed) 
(/!> along path B. The dashed line corresponds to the result at 
(l> — (j)c — 0.5265. Because (jtimilj ^) does not vary much with 
(j) we have not labelled the various curves with <j) 
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FIG. 13. </>"^(g ~ 4.7, w) at = 0.53, T* = 0.04 for 7 = ^ 
(solid), 1 (dotted) and 10 (dashed) 




2 4 6 8 10 12 14 16 18 20 
0) 



FIG. 14. g-dcpcndence of x"m('i')'^) at (/> = 0.381, 
T* = 0.04 for q ~ 2.0 (solid), q ~ 4.7 (dotted), q ~ 5.3 
(dashed), q ~ 6.5 (long dashed) and q ~ 9.8 (dot dashed) 
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FIG. 16. a - dependence of x'oo{<lmax,i^) at 4> = 0.381, 
T* = 0.04 for a = 0.06 (dashed), 0.065 (solid) and 0.068 
(dotted) 
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FIG. 15. g-dependence of (g, a;) at <;A = 0.53, T* = 0.04 
for q ~ 2.0 (solid), q ~ 4.7 (dotted), q ~ 5.3 (dashed), g ~ 6.5 
(long dashed) and q ~ 9.8 (dot dashed) 
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FIG. 17. (a) Temperature dependence of the position uiop 
(o) and height hop (□) (multiplied by 10) of the orientational 
peak for q ~ 4.7 and (j}D = 0.525 along path D. (b) De- 
pendence on the packing fraction of the position ujop (o) and 
height hop (□) (multiplied by 10) of the orientational peak for 
q ~ 4.7 and TX = 0.04 along path A 
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FIG. 18. Phase diagram for the glass transition between 
the phases I, II and III for the Bosse-Krieger model (p9|), (60) 
and (|^) with xo = 0.15. A' and B' denote paths along which 
we have investigated the dynamics 
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FIG. 19. 00 (t) for the Bosse-Krieger model for T = 0.7, 
4> = 3.9995, X = 0.15, SIq = 1, i^o = 10, Vi = 1, and different 
Q,x = 10""/^ n = 0,1,2,3,4 from left to right. jP denotes 
the critical non-ergodicity parameter of the _F2-model 
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CO 

FIG. 20. Xa(w) for fio = 1, f^i = 1/10, i^o = /^i = 0, 
X = 0.15 and T = 0.2 along path A' for = (1 + e)(?!)c, 
</>c ~ 2.8931318, £ = ±e"", n = 1, 2, 3 (sohd lines) and (i> = (t>c 
(dashed line). Because xliyi) does not vary much with we 
have not labelled the various curves with n 
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FIG. 21. x"(t^) for Ho = 1, fii = 1/10, va = = 0, 
X = 0.15 and T = 0.6 along path B' for = (l + e)(?!)c, 0c = 4, 
£ = -e"" for £ < and e = e""/^ for £ > 0, n = 1, 2, 3 (solid 
lines) and 4> = 4'c (dashed line). Because Xii'^) does not vary 
much with we have not labelled the various curves with n 
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FIG. 22. <;/>;,'(aj) for = 1, Qi = 1/10, I'd = i^i = 
a; = 0.15 and T = 0.2 along path A' for (j> = {1 + £)q 
cj>c ~ 2.8931318, e = ±e"", n = 1, 2, 3 (solid lines) and = 
(dashed line) 



FIG. 23. (t>a{i^) for fio = 1, tti 
a; = 0.15 and T = 0.6 along path B' for 
£ — — e^" for £ < and e = e~ 
lines) and ({> = <t>c (dashed line) 



1/10, uo = ui = 0, 
b= {l + e)(l)c, (jic = 4, 
for £ > 0, n = 1,2,3 (solid 
Because 4>i{'-^) does not vary 



much with A we have not labelled the various curves with n 
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FIG. 24. x'oi^) (jjoi^) for a reduced Bosse-Krieger 

model Qo = 1, fli = 1/10, uo = 0, ui = 0.19, x = 0.15 
and T = 0.6 (?!) = (1 + £)(l>c, (t)c = 4,, £ = -e~", for e < 



and £ 



-n/2 



for e > 0, n = 1,2,3 (solid lines) and 



(dashed line) 
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